Eccentricity fluctuations from the Color Glass Condensate at RHIC and LHC 
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In this brief note, we determine the fluctuations of the initial eccentricity in heavy-ion collisions 
caused by fluctuations of the nucleon configurations. This is done via a Monte-Carlo implementation 
of a Color Glass Condensate fct-factorization approach. The eccentricity fluctuations are found to 
nearly saturate elliptic flow fluctuations measured recently at RHIC. Extrapolations to LHC energies 
are shown. 
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PACS numbers: 12.38.Mh,24.85.+p,25.75.Ld,25.75.-q 

I. INTRODUCTION 

In a high energy non-central heavy ion collision the 
asymmetry of the coordinate space, the overlap area, is 
transferred into an asymmetry in momentum space, and 
measured as the elliptic flow V2 = (cos(20)). The initial 
asymmetry in coordinate space is characterized by the 
eccentricity, 



(1) 



where the brackets (...) indicate an average over the 
transverse plane, using some appropriate weight. Here, 
we use the number density of produced gluons. 

In ideal hydrodynamics, assuming a short thermaliza- 
tion time, the final elliptic flow is proportional to the 
initial eccentricity V2 — c e. The proportionality con- 
stant depends on the equation of state but is roughly 
c = 0.2 % 

Fluctuations of the eccentricity therefore should trans- 
late into fluctuations of the elliptic flow Q. Recently, 
these V2 fluctuations have been measured by the PHO- 
BOS and the STAR collaborations 0,(3. 

In this note, we examine the fluctuations of e based 
on the Monte Carlo KLN model introduced in Ref. [f| 
and compare to standard Glauber-model results (see, for 
example @, 0])- 



II. IMPROVEMENTS IN THE MC-KLN MODEL 

In Ref. [|| we introduced a Monte Carlo implemen- 
tation of the Kharzeev-Levin-Nardi (MC-KLN) [1| ap- 
proach to particle production in heavy ion collisions. 
Gluon production is calculated individually for each con- 
figuration of nucleons in the colliding nuclei. Thanks 
to the implementation of perturbative gluon saturation 
in this approach, the multiplicity can be determined via 
the well-known fct-factorization formula Q without the 
need to introduce infrared cutoffs (and additional mod- 
els for the soft regime). The saturation scale is taken to 
be proportional to the local density of nucleons which, 
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FIG. 1: (Color online) Multiplicity for Au+Au and Cu+Cu 
collisions at 200 GeV and PbPb collisions at 5500 GeV. The 
data is from the PHOBOS collaboration [|, L3- 



in turn, is measured by counting nucleons in a given 
sampling area. However, if the radius of the sampling 
area is r max = y <J- m ei /ft, one overestimates the interac- 
tion probability especially in the periphery, since nucleon 
pairs can have a distance up to 2r max . Therefore, we im- 
proved on our previous model by rejecting those pairs 
with r > r max . In the p + p limit this results in an addi- 
tional factor 0.58 which is very close to the value found 
in Refs. 0, HJ by accounting for the difference between 
the inelastic and the geometric cross section of a nucleon. 
We further assume here that a- me i = 42 mb at full RHIC 
energy (^/snn = 200 GeV), and a ine i = 66 mb at LHC 
energy (^snn = 5500 GeV). 

This refined treatment allows for an excellent descrip- 
tion of the charged multiplicity at RHIC over the en- 
tire range of centralities (for both Cu and Au nuclei), 
essentially down to p + p collisions. Fig. [1] depicts our 
results for full RHIC energy, as well as an extrapola- 
tion to Pb+Pb collisions at LHC energy. Since there is 
some uncertainty regarding the evolution of the satura- 
tion scale, we show results for both fixed coupling evo- 
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FIG. 2: (Color online) Relative fluctuations of the eccentricity 
as a function of centrality in Au+Au / Pb+Pb collisions. 



lution, Q 2 S = Q^ (xo/x) x with A = 0.28, and running 
coupling evolution of Q 2 S (see e.g. Q). For the latter 
case, the initial condition Q s ,o and xq was set such that 
at RHIC energy Q s agrees with previous estimates. 



account for the growth of dN/dr]/N part with centrality 
seen in Fig. [1]). These measurements are rather difficult, 
and therefore the error bars are quite large, as is the 
discrepancy between experiments, especially at high cen- 
tralizes where neither the Glauber model, nor the CGC 
result can be ruled out. For semi-central collisions, the 
CGC predicts somewhat lower relative fluctuations than 
the Glauber model. We note that a £ art itself is quite in- 
dependent of the underlying model and energy. The main 
reason for the lower relative eccentricity fluctuations in 
the MC-KLN model is the larger average eccentricity for 
semi-central Au+Au collisions in this approach, see the 
discussion in refs. [EG3- 

To check for other possible sources of fluctuations in 
the participant eccentricity, we implemented additional 
Poissonian (uncorrclated) fluctuations of the number of 
gluons produced at a given point in the transverse plane. 
These may arise, for example, from fluctuations of the 
gluon evolution ladders. However, we found that they did 
not noticeably affect er £part . One should also keep in mind 
that so-called non-flow effects may increase fluctuations 
of the measured vi- Moreover, hydrodynamic fluctua- 
tions may contribute to well Hence, a V2 /v 2 
should be viewed only as an upper limit for <7 epart /£part- 



III. FLUCTUATIONS OF THE INITIAL 
ECCENTRICITY 

The participant eccentricity e par t, which corrects for 
fluctuations of the major axes and of the center of mass 
of the overlap region, is defined by: 

v /( f Tg- t Tg) 2 +4< 

The fluctuations of this variable for a given centrality 
class (here defined by the number of participants) are 
determined via 

= \J (Epart) - (Epart) 2 ■ (3) 

Fig.JE] shows the result together with data from PHO- 
BOS [J and STAR [U, and a simple Glauber model, 
where the number density of gluons scales with the num- 
ber of participants N part (note that this model fails to 



IV. SUMMARY 

We have calculated the fluctuations of the initial ec- 
centricity within a simple Glauber model and for a Color 
Glass Condensate approach which includes fluctuations 
in the positions of the hard sources (nuclcons). Both 
models predict eccentricity fluctuations which nearly sat- 
urate the experimentally measured fluctuations of the el- 
liptic flow. The CGC approach gives slightly lower rela- 
tive fluctuations than the Glauber model, which is largely 
due to a higher average eccentricity £ par t- Their magni- 
tude at LHC energy is similar. 
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